Identification of genotype-biochemical phenotype correlations associated with fructose 1,6-bisphosphatase deficiency.

Fructose-1,6-bisphosphatase (FBPase) deficiency, caused by an FBP1 mutation, is an autosomal recessive disorder characterized by hypoglycemic lactic acidosis. Due to the rarity of FBPase deficiency, the mechanism by which the mutations cause enzyme activity loss still remains unclear. Here we identify compound heterozygous missense mutations of FBP1, c.491G>A (p.G164D) and c.581T>C (p.F194S), in an adult patient with hypoglycemic lactic acidosis. The G164D and F194S FBP1 mutants exhibit decreased FBP1 protein expression and a loss of FBPase enzyme activity. The biochemical phenotypes of all previously reported FBP1 missense mutations in addition to G164D and F194S are classified into three functional categories. Type 1 mutations are located at pivotal residues in enzyme activity motifs and have no effects on protein expression. Type 2 mutations structurally cluster around the substrate binding pocket and are associated with decreased protein expression due to protein misfolding. Type 3 mutations are likely nonpathogenic. These findings demonstrate a key role of protein misfolding in mediating the pathogenesis of FBPase deficiency, particularly for Type 2 mutations. This study provides important insights that certain patients with Type 2 mutations may respond to chaperone molecules.

Carcinogenic potential in regenerated mucosa after endoscopic resection of esophageal squamous cell carcinoma.

BACKGROUND AND AIM: Little is known about genetic mutations in the regenerated mucosa (RM) after endoscopic resection (ER) of esophageal carcinoma. Thus, this study investigates the status of genetic variation in RM after ER of esophageal squamous cell carcinoma (ESCC). METHODS: The study cohort included 19 patients with ESCC. We used an esophageal carcinoma panel to identify target sequences for squamous cell carcinoma (SCC), background mucosa (BM), and RM after ER of ESCC. We used OncoKB to check whether each mutation was a putative driver. RESULTS: We identified 77 mutations of 32 genes in SCC, 133 mutations of 34 genes in BM, and 100 mutations of 29 genes in RM. Putative driver mutations were identified in 20 mutations in 14 cases in SCC, 16 mutations in 10 cases in BM, and 7 mutations in 11 cases in RM. The rate of putative driver mutations to total mutations was significantly lower in RM (26% in SCC vs 12% in BM vs 7% in RM, P = 0.009). Additionally, the rate of cases with TP53 putative driver mutations was significantly lower in RM (63% in SCC vs 37% in BM vs 16% in RM, P = 0.011). The percentage of putative driver mutations and the percentage of cases with a putative driver of TP53 were significantly lower in RM. CONCLUSION: Esophageal RM after ER of ESCC could have a lower risk of carcinogenesis.

Extension of bacterial rDNA sequencing for simultaneous methylation detection and its application in microflora analysis.

Although polymerase chain reaction (PCR) amplification and sequencing of the bacterial 16S rDNA region has numerous scientific applications, it does not provide DNA methylation information. Herein, we propose a simple extension for bisulfite sequencing to investigate 5-methylcytosine residues in the bacterial 16S rDNA region from clinical isolates or flora. Multiple displacement amplification without DNA denaturation was used to preferentially pre-amplify single-stranded bacterial DNA after bisulfite conversion. Following the pre-amplification, the 16S rDNA region was analyzed using nested bisulfite PCR and sequencing, enabling the simultaneous identification of DNA methylation status and sequence data. We used this approach (termed sm16S rDNA PCR/sequencing) to identify novel methylation sites and a methyltransferase (M. MmnI) in Morganella morganii and different methylation motifs among Enterococcus faecalis strains from small volumes of clinical specimens. Further, our analysis suggested that M. MmnI may be correlated to erythromycin resistance. Thus, sm16S rDNA PCR/sequencing is a useful extension method for analyzing the DNA methylation of 16S rDNA regions in a microflora, providing additional information not provided by conventional PCR. Given the relationship between DNA methylation status and drug resistance in bacteria, we believe this technique can be effectively applied in clinical sample testing.

Validity of pathological diagnosis for early colorectal cancer in genetic background.

BACKGROUND: This study aimed to investigate the validity of pathological diagnosis of early CRC (E-CRC) from the genetic background by comparing data of E-CRC to colorectal adenoma (CRA) and The Cancer Genome Atlas (TCGA) on advanced CRC (AD-CRC). METHODS: TCGA data on AD-CRC were studied in silico, whereas by next-generation sequencer, DNA target sequences were performed for endoscopically obtained CRA and E-CRC samples. Immunohistochemical staining of mismatch repair genes and methylation of MLH1 was also performed. The presence of oncogenic mutation according to OncoKB for the genes of the Wnt, MAPK, and cell-cycle-signaling pathways was compared among CRA, E-CRC, and AD-CRC. RESULTS: The study included 22 CRA and 30 E-CRC lesions from the Chiba University Hospital and 212 AD-CRC lesions from TCGA data. Regarding the number of lesions with driver mutations in the Wnt and cell-cycle-signaling pathways, E-CRC was comparable to AD-CRC, but was significantly greater than CRA. CRA had significantly more lesions with a driver mutation for the Wnt signaling pathway only, versus E-CRC. CONCLUSIONS: In conclusion, the definition of E-CRC according to the Japanese criteria had a different genetic profile from CRA and was more similar to AD-CRC. Based on the main pathway, it seemed reasonable to classify E-CRC as adenocarcinoma. The pathological diagnosis of E-CRC according to Japanese definition seemed to be valid from a genetic point of view.

MicroRNA-874 targets phosphomevalonate kinase and inhibits cancer cell growth via the mevalonate pathway.

The microRNA (miR) miR-874, a potential tumour suppressor, causes cell death via target gene suppression in various cancer types. Mevalonate pathway inhibition also causes cell death in breast cancer. However, the relationship between the mevalonate pathway and miR-874-induced apoptosis or its association with the tumour suppressor p53 has not been elucidated. We identified phosphomevalonate kinase (PMVK), a key mevalonate pathway enzyme, and sterol regulatory element-binding factor 2 (SREBF2), the master cholesterol biosynthesis regulator, as direct miR­874 targets. Next-generation sequencing analysis revealed a significant miR-874-mediated downregulation of PMVK and SREBF2 gene expression and p53 pathway enrichment. Luciferase reporter assays showed that miR-874 directly regulated PMVK and SREBF2. miR-874-induced apoptosis was p53 dependent, and single-cell RNA sequencing analysis demonstrated that miR-874 transfection resulted in apoptosis and p53 pathway activation. Downregulation of PMVK expression also caused cell cycle arrest and p53 pathway activation, which was rescued by geranylgeranyl pyrophosphate (GGPP) supplementation. Analysis of The Cancer Genome Atlas (TCGA) database indicated a negative correlation between miR-874 and PMVK expression and between miR-874 and SREBF2 expression. These findings suggest that miR-874 suppresses the mevalonate pathway by targeting SREBF2 and PMVK, resulting in GGPP depletion, which activates the p53 pathway and promotes cycle arrest or apoptosis.

Proteogenomic landscape and clinical characterization of GH-producing pituitary adenomas/somatotroph pituitary neuroendocrine tumors.

The clinical characteristics of growth hormone (GH)-producing pituitary adenomas/somatotroph pituitary neuroendocrine tumors (GHomas/somatotroph PitNETs) vary across patients. In this study, we aimed to integrate the genetic alterations, protein expression profiles, transcriptomes, and clinical characteristics of GHomas/somatotroph PitNETs to identify molecules associated with acromegaly characteristics. Targeted capture sequencing and copy number analysis of 36 genes and nontargeted proteomics analysis were performed on fresh-frozen samples from 121 sporadic GHomas/somatotroph PitNETs. Targeted capture sequencing revealed GNAS as the only driver gene, as previously reported. Classification by consensus clustering using both RNA sequencing and proteomics revealed many similarities between the proteome and the transcriptome. Gene ontology analysis was performed for differentially expressed proteins between wild-type and mutant GNAS samples identified by nontargeted proteomics and involved in G protein-coupled receptor (GPCR) pathways. The results suggested that GNAS mutations impact endocrinological features in acromegaly through GPCR pathway induction. ATP2A2 and ARID5B correlated with the GH change rate in the octreotide loading test, and WWC3, SERINC1, and ZFAND3 correlated with the tumor volume change rate after somatostatin analog treatment. These results identified a biological connection between GNAS mutations and the clinical and biochemical characteristics of acromegaly, revealing molecules associated with acromegaly that may affect medical treatment efficacy.

Development of multiplex S-gene-targeted RT-PCR for rapid identification of SARS-CoV-2 variants by extended S-gene target failure.

BACKGROUND: Tracking SARS-CoV-2 variants of concern (VOC) by genomic sequencing is time-consuming. The rapid screening of VOCs is necessary for clinical laboratories. In this study, we developed a rapid screening method based on multiplex RT-PCR by extended S-gene target failure (eSGTF), a false negative result caused by S-gene mutations. METHODS: Three S-gene target (SGT) regions (SGT1, codons 65-72; SGT2, codons 152-159; and SGT3, codons 370-377) and an N-gene region (for internal control) were detected in single-tube. Four types of VOC (Alpha, Delta, Omicron BA.1, and Omicron BA.2) are classified by positive/negative patterns of 3 S-gene regions (eSGTF pattern). RESULTS: The eSGTF patterns of VOCs were as follows (SGT1, SGT2, SGT3; P, positive; N, negative): Alpha, NPP; Delta, PNP; Omicron BA.1, NPN pattern; and Omicron BA.2, PPN. As compared with the S-gene sequencing, eSGTF patterns were identical to the specific VOCs (concordance rate = 96.7%, N = 206/213). Seven samples with discordant results had a minor mutation in the probe binding region. The epidemics of VOCs estimated by eSGTF patterns were similar to those in Japan. CONCLUSIONS: Multiplex RT-PCR and eSGTF patterns enable high-throughput screening of VOCs. It will be useful for the rapid determination of VOCs in clinical laboratories.

DRD2 Taq1A Polymorphism-Related Brain Volume Changes in Parkinson's Disease: Voxel-Based Morphometry.

Taq1A polymorphism is a DRD2 gene variant located in an exon of the ANKK1 gene and has an important role in the brain's dopaminergic functions. Some studies have indicated that A1 carriers have an increased risk of developing Parkinson's disease (PD) and show poorer clinical performance than A2 homo carriers. Previous studies have suggested that A1 carriers had fewer dopamine D2 receptors in the caudate and increased cortical activity as a compensatory mechanism. However, there is little information about morphological changes associated with this polymorphism in patients with PD. The study's aim was to investigate the relationship between brain volume and Taq1A polymorphism in PD using voxel-based morphometry (VBM). Based on Taq1A polymorphism, 103 patients with PD were divided into two groups: A1 carriers (A1/A1 and A1/A2) and A2 homo carriers (A2/A2). The volume of the left prefrontal cortex (PFC) was significantly decreased in A2 homo carriers compared to A1 carriers. This finding supports the association between Taq1A polymorphism and brain volume in PD and may explain the compensation of cortical function in A1 carriers with PD.

Tracking SARS-CoV-2 variants by entire S-gene analysis using long-range RT-PCR and Sanger sequencing.

INTRODUCTION: Genomic surveillance of the SARS-CoV-2 virus is important to assess transmissibility, disease severity, and vaccine effectiveness. The SARS-CoV-2 genome consists of approximately 30 kb single-stranded RNA that is too large to analyze the whole genome by Sanger sequencing. Thus, in this study, we performed Sanger sequencing following long-range RT-PCR of the entire SARS-CoV-2 S-gene and analyzed the mutational dynamics. METHODS: The 4 kb region, including the S-gene, was amplified by two-step long-range RT-PCR. Then, the entire S-gene sequence was determined by Sanger sequencing. The amino acid mutations were identified as compared with the reference SARS-CoV-2 genome. RESULTS: The S:D614G mutation was found in all samples. The R.1 variants were detected after January 2021. Alpha variants started to emerge in April 2021. Delta variants replaced Alpha in July 2021. Then, Omicron variants were detected after December 2021. These mutational dynamics in samples collected in the Chiba University Hospital were similar to those in Japan. CONCLUSION: The emergence of variants of concern (VOC) has been reported by the entire S-gene analysis. As the VOCs have unique mutational patterns of the S-gene region, analysis of the entire S-gene will be useful for molecular surveillance of the SARS-CoV-2 in clinical laboratories.

Genetic profiles of Barrett's esophagus and esophageal adenocarcinoma in Japanese patients.

The genetic characteristics of Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC) in the Japanese population is unclear. This study aims to investigate the genetic characteristics from nondysplastic BE (NDBE) to early EAC in Japan. Clinical information was collected. Moreover, the genetic profile of NDBE without concurrent dysplasia, early EAC, and surrounding BE were also investigated using endoscopic biopsy samples and formalin-fixed, paraffin-embedded specimens from Japanese patients by targeted next-generation sequencing. Immunohistochemical staining for p53 was also performed for EAC lesions. Targeted NGS was performed for 33 cases with 77 specimens. No significant difference exists in the NDBE group between the number of putative drivers per lesion in the short-segment Barrett's esophagus (SSBE) and long-segment Barrett's esophagus (LSBE) [0 (range, 0-1) vs. 0 (range, 0-1). p = 1.00]. TP53 putative drivers were found in two patients (16.7%) with nondysplastic SSBE. TP53 was the majority of putative drivers in both BE adjacent to EAC and EAC, accounting for 66.7% and 66.7%, respectively. More putative drivers per lesion were found in the EAC than in the NDBE group [1 (range, 0-3) vs. 0 (range, 0-1). p < 0.01]. The genetic variants of TP53 in the Japanese early EAC were similar to those in western countries. However, TP53 putative drivers were detected even in Japanese patients with nondysplastic SSBE. This is significant because such nondysplastic SSBE might have higher risk of progressing to high-grade dysplasia or EAC. The risks of progression may not be underestimated and appropriate follow-ups may be necessary even in patients with SSBE.Trial registration: This study was registered at the University Hospital Medical Information Network (UMIN000034247).

Comprehensive mutational analysis of background mucosa in patients with Lugol-voiding lesions.

Somatic mutations including the background mucosa in patients with Lugol-voiding lesions (LVLs) are still not well known. The aim of this study was to evaluate the somatic mutations of the background mucosa in patients with LVLs (Squamous cell carcinoma (SCC), intraepithelial neoplasia (IN), and hyperplasia). Twenty-five patients with LVLs (9 with SCC, 6 with IN, and 10 with hyperplasia) were included. A targeted sequence was performed for LVLs and background mucosa using an esophageal cancer panel. Each mutation was checked whether it was oncogenic or not concerning OncoKB. In LVLs, TP53 was the most dominant mutation (80%). Furthermore, 72% of TP53 mutations was putative drivers. In background mucosa, NOTCH1 was the most dominant mutation (88%) and TP53 was the second most dominant mutation (48%). Furthermore, 73% of TP53 mutations and 8% of NOTCH1 mutations were putative drivers. Putative driver mutations of TP53 had significantly higher allele frequency (AF) in SCC than in IN and hyperplasia. Conversely, putative driver mutations of NOTCH1 did not have a significant accumulation of AF in the progression of carcinogenesis. Furthermore, in SCC, AF of TP53 mutations was significantly higher in LVLs than in background mucosa, but not in IN and hyperplasia. Regarding NOTCH1, a significant difference was not observed between LVLs and background mucosa in each group. The background mucosa in patients with LVLs already had putative driver mutations such as TP53 and NOTCH1. Of these two genes, TP53 mutation could be the main target gene of carcinogenesis in esophageal SCC. Clinical Trials registry: UMIN000034247.

TAS4464, a NEDD8-activating enzyme inhibitor, activates both intrinsic and extrinsic apoptotic pathways via c-Myc-mediated regulation in acute myeloid leukemia.

TAS4464, a potent, selective small molecule NEDD8-activating enzyme (NAE) inhibitor, leads to inactivation of cullin-RING E3 ubiquitin ligases (CRLs) and consequent accumulations of its substrate proteins. Here, we investigated the antitumor properties and action mechanism of TAS4464 in acute myeloid leukemia (AML). TAS4464 induced apoptotic cell death in various AML cell lines. TAS4464 treatments resulted in the activation of both the caspase-9-mediated intrinsic apoptotic pathway and caspase-8-mediated extrinsic apoptotic pathway in AML cells; combined treatment with inhibitors of these caspases markedly diminished TAS4464-induced apoptosis. In each apoptotic pathway, TAS4464 induced the mRNA transcription of the intrinsic proapoptotic factor NOXA and decreased that of the extrinsic antiapoptotic factor c-FLIP. RNA-sequencing analysis showed that the signaling pathway of the CRL substrate c-Myc was enriched after TAS4464 treatment. Chromatin immunoprecipitation (ChIP) assay revealed that TAS4464-induced c-Myc bound to the PMAIP1 (encoding NOXA) and CFLAR (encoding c-FLIP) promoter regions, and siRNA-mediated c-Myc knockdown neutralized both TAS4464-mediated NOXA induction and c-FLIP downregulation. TAS4464 activated both caspase-8 and caspase-9 along with an increase in NOXA and a decrease in c-FLIP, resulting in complete tumor remission in a human AML xenograft model. These findings suggest that NAE inhibition leads to anti-AML activity via a novel c-Myc-dependent apoptosis induction mechanism.

Comprehensive Analysis of Barrett's Esophagus: Focused on Carcinogenic Potential for Barrett's Cancer in Japanese Patients.

BACKGROUND/AIM: Barrett's esophagus (BE) is a precursor of esophageal adenocarcinoma (EAC). Therefore, an accurate diagnosis of BE is important for the subsequent follow-up and early detection of EAC. However, the definitions of BE have not been standardized worldwide; columnar-lined epithelium (CLE) without intestinal metaplasia (IM) and/or < 1 cm is not diagnosed as BE in most countries. This study aimed to clarify the malignant potential of CLE without IM and/or < 1 cm genetically. METHOD: A total of 96 consecutive patients (including nine patients with EAC) who had CLE were examined. Biopsies for CLE were conducted, and patients were divided into those with IM and > 1 cm (Group A) and those without IM and/or < 1 cm (Group B). Malignant potential was assessed using immunochemical staining for p53. Moreover, causative genes were examined using next-generation sequencing (NGS) on ten patients without Helicobacter pylori infection and without atrophic gastritis. RESULT: Of the 96 patients, 66 were in Group B. The proportion of carcinoma/dysplasia in Group A was significantly higher than that in Group B (26.7% in Group A and 1.5% in Group B; p < 0.01). However, one EAC patient was found in Group B. In the immunostaining study for non-EAC patients, an abnormal expression of p53 was not observed in Group A, whereas p53 loss was observed in three patients (4.6%) in Group B. In the NGS study, a TP53 mutation was found in Group B. CONCLUSION: CLE without IM and/or < 1 cm has malignant potential. This result suggests that patients with CLE as well as BE need follow-up.

Monoamine oxidase B rs1799836 G allele polymorphism is a risk factor for early development of levodopa-induced dyskinesia in Parkinson's disease.

BACKGROUND: Dopamine replacement therapy is an established treatment for motor symptoms of Parkinson's disease, but its long-term use is often limited by the eventual development of motor complications, including levodopa-induced dyskinesia. Genetic background, particularly polymorphisms of dopamine metabolism genes, may affect the occurrence of dyskinesia in Parkinson's disease patients. METHODS: We investigated polymorphisms of dopamine metabolism genes, including catechol-O-methyltransferase, monoamine oxidase B, dopamine beta-hydroxylasedopamine, dopamine receptors D1, D2, and D3, and dopamine transporter, in 110 patients with Parkinson's disease. Cox proportional hazards regression was used to detect associations between genotypes and levodopa-induced dyskinesia. RESULTS: Monoamine oxidase B rs1799836 was the only polymorphism correlated with risk of dyskinesia. Patients with an AG or GG genotype were more likely to have dyskinesia than those with an AA genotype (adjusted hazard ratio, 3.41; 95% confidence interval, 1.28-9.10). Also, Kaplan-Meier curves demonstrated that patients with an AG or GG genotype developed dyskinesia earlier than those with an AA genotype (log-rank test, p = .004). CONCLUSIONS: In Parkinson's disease patients, the monoamine oxidase B rs1799836 G allele is associated with a greater likelihood of developing dyskinesia than the A allele, possibly due to its association with lower monoamine oxidase B activity in the brain. Thus, detection of monoamine oxidase B polymorphisms may be useful for determining the optimal dosing of antiparkinson medications.

Highly sensitive detection of SARS-CoV-2 RNA by multiplex rRT-PCR for molecular diagnosis of COVID-19 by clinical laboratories.

BACKGROUND: The detection of SARS-CoV-2 RNA by real-time reverse transcription-polymerase chain reaction (rRT-PCR) is used to confirm the clinical diagnosis of COVID-19 by molecular diagnostic laboratories. We developed a multiplex rRT-PCR methodology for the detection of SARS-CoV-2 RNA. METHODS: Three genes were used for multiplex rRT-PCR: the Sarbecovirus specific E gene, the SARS-CoV-2 specific N gene, and the human ABL1 gene as an internal control. RESULTS: Good correlation of Cq values was observed between the simplex and multiplex rRT-PCR methodologies. Low copies (<25 copies/reaction) of SARS-CoV-2 RNA were detected by the novel multiplex rRT-PCR method. CONCLUSION: The proposed multiplex rRT-PCR methodology will enable highly sensitive detection of SARS-CoV-2 RNA, reducing reagent use and cost, and time required by clinical laboratory technicians.

Central blood pressure and pulse wave velocity in young and middle-aged Japanese adults with isolated systolic hypertension.

Isolated systolic hypertension (ISH), defined as systolic blood pressure (SBP) ≥140 mmHg and diastolic BP (DBP) <90 mmHg, is a common type of hypertension among young men. This study aimed to investigate the clinical characteristics, central blood pressure, and arterial stiffness of young and middle-aged Japanese individuals with ISH. A total of 432 male participants, aged 18-49 years, were classified into six subgroups: optimal BP (SBP <120 mmHg and DBP <80 mmHg), high-normal BP (SBP 120-129 mmHg and DBP <80 mmHg), high-BP (SBP 130-139 mmHg and/or DBP 80-89 mmHg), ISH (SBP ≥140 mmHg and DBP <90 mmHg), isolated diastolic hypertension (IDH) (SBP <140 mmHg and DBP ≥90 mmHg), and systolic and diastolic hypertension (SDH) (SBP ≥140 mmHg and DBP ≥90 mmHg). Participants with ISH had a greater body mass index (BMI) and waist circumference than the optimal BP participants but were more likely to be physically active than the IDH and SDH participants. The central SBP of the ISH subgroup was higher than that of the optimal/high-normal/high-BP subgroups and lower than that of the SDH subgroup. The carotid-femoral pulse wave velocity (cfPWV) of the ISH subgroup was higher than that of the optimal and high-normal BP subgroups and lower than that of the SDH subgroup after adjusting for age, heart rate, BMI, and physical activity. These differences disappeared after further adjustment for central mean arterial pressure. In conclusion, the central SBP of Japanese men with ISH was greater than that of Japanese men with optimal/high-normal/high-BP, but the progression of arterial stiffness was unclear.

A simple method for sequencing the whole human mitochondrial genome directly from samples and its application to genetic testing.

Next-generation sequencing (NGS) is a revolutionary sequencing technology for analyzing genomes. However, preprocessing methods for mitochondrial DNA (mtDNA) sequencing remain complex, and it is required to develop an authenticated preprocessing method. Here, we developed a simple and easy preprocessing method based on isothermal rolling circle mtDNA amplification using commercially available reagents. Isothermal amplification of mtDNA was successfully performed using both nanoliter quantities of plasma directly and 25 ng of total DNA extracted from blood or tissue samples. Prior to mtDNA amplification, it was necessary to treat the extracted total DNA with Exonuclease V, but it was not required to treat plasma. The NGS libraries generated from the amplified mtDNA provided sequencing coverage of the entire human mitochondrial genome. Furthermore, the sequencing results successfully detected heteroplasmy in patient samples, with called mutations and variants matching those from previous, independent, Sanger sequencing analysis. Additionally, a novel single nucleotide variant was detected in a healthy volunteer. The successful analysis of mtDNA using very small samples from patients is likely to be valuable in clinical medicine, as it could reduce patient discomfort by reducing sampling-associated damage to tissues. Overall, the simple and convenient preprocessing method described herein may facilitate the future development of NGS-based clinical and forensic mtDNA tests.

Multiplex PCR and multicolor probes melting for the simultaneous detection of five UGT1A1 variants.

The multiplex PCR melting analysis method was developed for detecting the five UGT1A1 variants. Multiplexing was achieved using color probes and Tm. The probes for *28/*6, *27, *29, and *7 were discriminated by colors. Although the probes for *28 and *6 had the same colors, their variants were clearly discriminated by probe Tm. The allelic frequencies of each genotype were 0.12 for *28, 0.19 for *6, 0.02 for *27, 0.0 for *29, and 0.005 for *7. We developed a multiplex PCR melting analysis method, which will be useful in molecular diagnostics and pharmacogenetic analyses in clinical laboratories.

Evaluation of analytical factors associated with targeted MEFV gene sequencing using long-range PCR/massively parallel sequencing of whole blood DNA for molecular diagnosis of Familial Mediterranean fever.

BACKGROUND: Long-range PCR (LR-PCR) is used to enrich the target regions of the genome. This study aimed to establish the pipeline of targeted gene sequencing using LR-PCR and massively parallel sequencing (MPS). METHODS: The 14-kb-long MEFV gene, including the entire coding exons, was selected as a target gene and amplified using LR-PCR. The evaluated analytical factors were as follows: LR-PCR conditions, three types of post-PCR cleanup methods, and two types of MPS library preparation methods. RESULTS: With regard to LR-PCR conditions, Tks Gflex DNA polymerase at 7-min (30-s/kb) annealing/extension with 100-ng genomic DNA input had the highest yield. Regarding post-PCR purification methods, the magnetic beads-based method had high recovery and purity. In the MPS library preparation methods, the ligation-based method had a higher base coverage in the target (94.58%), uniformity of base coverage (99.95%), and target bases with no strand bias (97.40%). The exonic variants determined by Sanger sequencing were detected by both ligation- and transposon-based methods. CONCLUSIONS: Various analytical factors were evaluated, and the pipeline of targeted gene sequencing using LR-PCR and MPS was established. These data can enable the optimization of targeted gene sequencing using LR-PCR and MPS in the clinical laboratory.